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Abstract 

The 2-point functions for AS = 1 current- current and QCD-penguin operators, as well 
as for the AS = 2 operator, are calculated at the next-to-leading order. The calculation 
is performed in two different renormalization schemes for 75, and the compatibility of the 
results obtained in the two schemes is verified. The scale- and scheme-invariant combina- 
tions of spectral 2-point functions and corresponding Wilson-coefficients are constructed 
and analyzed. For AS = 1, the QCD corrections to the CP-conserving part, dominated by 
current-current operators, are 40%-120% at = (1 — 3 GeV)^, whereas the correction to 
the imaginary part, mainly coming from the penguin operator Qq, are 100%-240%. The 
large size of the gluonic corrections to current-current operators provides a qualitative un- 
derstanding of the observed enhancement of A J = 1/2 transitions. In the AS = 2 sector 
the QCD corrections are quite moderate (~ —20%). 
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1 Introduction 



Recent years have witnessed considerable improvement in Standard Model calculations of 
non-leptonic weak decays. In particular, the effective Hamiltonians for flavour- changing 
AF = 1 as well as AF = 2 transitions were calculated at the next-to- leading 

order (NLO) in renormalization group (RG) improved perturbation theory. The neces- 
sary computations included the determination of 2-loop anomalous dimension matrices for 
current-current, QCD-penguin, and electroweak penguin operators [|I],0,|^||,y,|TO|,|Tl|. The 
benefit from such calculations is the following: 

• the size of the NLO short distance corrections to the coefficient functions was ob- 
tained, 

• this allowed an estimation of the scale down to which RG evolution is possible, before 
perturbation theory breaks down, to be around 1 GeV, 

• only a NLO calculation allows for a meaningful use of the QCD scale Aj^, extracted 
for example in deep inelastic scattering, r decay or LEP experiments, 

• and it made possible to attack the question of the dependence of coefficient functions 
on the renormalization scheme, e.g. the definition of 75 in an arbitrary space-time 
dimension, which first appears at the next-to-leading order. 

Sadly enough, this is not the whole story. In order to fully calculate the decay amplitude 
for a certain process, we also need to know the matrix elements of the operators appearing 
in the effective Hamiltonian, between the hadronic initial and final states. This part is 
much more difficult, for it involves non-perturbative dynamics at low energies. Methods to 
attempt this involved task include lattice gauge theory |12,1^, 14|, 1/A^-expansion |15,16 



|T^, chiral perturbation theory [jT8|,|T9|, QCD sum rules [^|2T|,^,^^,^,^, and mixed 



approaches involving functional integration of quark fields [^. A strategy to obtain the 
matrix elements for AS" = 1 decays at NLO as far as possible from experimental data was 
also advocated in ref. [^]. 

All these methods suffer from more or less severe drawbacks. Although the lattice 
could eventually become the ultimate tool for the calculation of matrix elements, the pre- 
cision of present lattice results is still very poor. 1/A^-expansion and chiral perturbation 
theory are not yet directly related to the fundamental QCD Lagrangian, and therefore, 
a sound matching of matrix elements calculated in one of these methods and the coeffi- 
cient functions, obtained in perturbation theory, is still not possible. The bridge between 
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the effective low-energy chiral Lagrangian and the underlying QCD theory can be built 



with functional bosonization techniques or using QCD sum rules to relate the two 



energy regimes [^,gT],^,g3|,|2^,g5|,|22l. Unfortunately, the present determinations are not 
very accurate. Finally, there is not enough experimental information to obtain all matrix 
elements, say for K —>■ tttt decays, from the phenomenological method of ref. |Q. 

The problem becomes much easier at the inclusive level, where the properties of the 
non-leptonic effective weak Hamiltonian can be analyzed within QCD [2^,^,^,^,^^ . 
Given, for instance, the short- distance AS* = 1 Hamiltonian 

= ^v^,v:,Y.c^{^^')Q^, (1.1) 

obtained through the operator product expansion, one considers the 2-point function 

This vacuum-to-vacuum correlator can be studied with perturbative QCD methods, al- 
lowing for a consistent combination of Wilson-coefficients Cj(/i^) and 2-point functions 
of the 4-quark operators, in such a way that the renormalization scheme and scale 
dependences exactly cancel (to the computed order). The associated spectral function 
^lm'^^^^^{q'^) is a quantity with definite physical information. It describes in an inclusive 
way how the weak Hamiltonian couples the vacuum to physical states of a given invariant 
mass. General properties like the observed enhancement of A/ =1/2 transitions can be 
then rigorously analyzed at the inclusive level. 

A detailed analysis of 2-point functions associated with AS* = 1 and AS* = 2 operators 
was presented in ref. p6[, where the 0{as) corrections to the corresponding correlators 
were calculated. The NLO corrections to the A/ = 1/2 2-point functions were found to 



be very large |]26|, confirming the QCD enhancement obtained in a previous approximate 



calculation [24]. The results of ref. |2^ were, however, incomplete because the NLO correc- 



tions to the Wilson-coefficients of penguin operators were still missing. With the progress 
achieved for the Wilson-coefficient functions mentioned above, we are now in a position to 
match matrix elements and coefficient functions consistently at NLO. 

To get a sensible result, we obviously need to use the same renormalization scheme 
conventions on both sides of the calculation. Unfortunately, for technical reasons, different 
bases of operators have been used in the 2-point function and Wilson-coefficient calcula- 
tions. In order to avoid ambiguities coming from the definition of 75 in 7^ 4 dimensions, a 
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set of colour-singlet 4-quark operators was used in ref. |]26| to perform the 2-point-function 
calculation; the computation was done with dimensional regularization and a naively anti- 
commuting 75 (NDR scheme). While that is fine in the leading logarithmic approximation, 
the basis of colour-singlet operators does not close under renormalization at the NLO. The 
Fierz-transformations, which are needed to relate some of the operators in the process of 
renormalization, are broken by 0{as) corrections, and additional contributions have to 
be taken into account. For this reason, we shall reconsider the calculation of the 2-point 
functions of 4-quark operators in this work. 

We shall use the same basis of operators which has been taken for the calculation of the 
Wilson-coefficients, so that we can directly incorporate the results of refs. [0,|^,H,|Tl|. The 
presence of colour-non-singlet operators in this basis gives rise to 75 complications in the 
2-point-function evaluation. Like for the calculation of the anomalous dimension matrices 
in refs. 0,^H,|Tl|, we shall perform the calculation in two different schemes for 75, to have 
explicit tests on our result. A direct computation with a naively anticommuting 75 is not 
possible, since two diagrams include traces of odd numbers of 75; but we shall show, how 
nevertheless a result in the NDR scheme can be obtained. For the second computation, 
the consistent definition of a non-anticommuting 75 in arbitrary dimensions according to 
't Hooft and Veltman (HV scheme) 0,^,0 is used. 

In sect. 2, we shall discuss the general structure of 2-point functions of 4-quark op- 
erators. As a first step towards the explicit calculation for the AS* = 1 case, in sect. 3, 
the 2-point functions of current-current operators are computed, and the full set including 
QCD-penguins is presented in sect. 4. A numerical analysis of the results is given in sect. 5. 
In sect. 6, we evaluate the 2-point function for the AS = 2 operators. A comparison with 



the results of ref. 26 , together with some concluding remarks, is finally given in sect. 7. 



2 General structure 

As a first step, let us discuss the general structure of the 2-point functions of 4-quark 
operators and their renormalization. The bare 2-point function \E'^(g^) is defined by 

^^(g2) = r y'(ixe*^^(0|r{g^(x)g^'(0)}|0). (2.1) 

Q{x) can either be a single operator, or a vector of 4-quark operators, in which case \E'^(g^) 
i= a symmetric matr ij. 

^For this general discussion, we shall assume '^^{q^) to be a matrix. 
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Keeping only relevant terms up to next-to-leading order in as, and working with dimen- 
sional regularization, the regularized, but yet unrenormalized 2-point function \l/^(g^,/i^) 
has an expansion in e {d = 4 + 2e), 

-q^\^' 1 r 



2\4 



(4vr) 



3e 



+ 



a. 



IT 



-Q^\^' 1 



4^2 



A + Be 



C + De + 



+ 0(ai 



(2.2) 



H is a renormalization scale in the MS scheme, that is, we have redefined the scale of 
dimensional regularization z/^ to be exp(7£;) fi'^/{47i), where je is Euler's constant, so that 
only poles in e have to be subtracted. Our main goal will be to calculate the four matrices 
A, B, C, and D. 

The renormalized 2-point function is given by 



-l\T 



(2.3) 



where Z is the renormalization matrix of the 4-quark operators, Q = Z Q, and Rjjs 
means that additional poles in e, stemming from the operator product, need to be sub- 
tracted. In the minimal subtraction scheme, Z has the general expansion 



zi^^ 



(2.4) 



k=l " n=l ^ 

The anomalous dimension matrix of 4-quark operators is defined through 
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7-1 dZ 
Z fi — 
djjL 



7 



(1) 



e=0 



^ + ^(2) f^f + ... . 
TT ^ -K ' 



Inserting the expansion for Z, eq. ( p.4|) , we find 
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(1) 



2Z\ 



(1) 



7 



(2) 



AZ\ 



(2) 



'(2) 



-Z^^\2Z\ 



(1) 



/9i 



(2.5) 



(2.6) 



Here, I3i = — {IIN — 2/)/6 is the leading coefficient of the /3-function, and / being 
the number of colours and fiavours respectively. For the operators which mediate AS* = 1 
transitions, and shall be our interest for most part of the paper, the leading order anomalous 



dimension matrix is known already since a long time P9|, p0i[31| ,|32|,p3|. However, the full 
next-to-leading order matrix has only been obtained recently [|l],0,|^,|^,|lT|. To make contact 
with the notation of refs. we note that 
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(1) 



-L (0) 

^ Ibjlw 



and 7(2) = — 7BJLW 



(2.7) 



4 



Using the results of eq. ( pl6| ) , together with eq. ( p^ ) , up to non-logarithmic corrections 
the renormalized 2-point function turns out to be 



where L = ln(— g //i ), and 



2\4 



X 



(Any 



D - 



AL + 



IT [2 



CL^ + XL 



(2.8) 
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Because of renormalizability, it also follows that 



C 



(2.9) 



(2.10) 



For the rest of this work, we shall only be concerned with the spectral function (the 
imaginary part of the 2-point function) which is directly related to physical quantities: 



$(s,yU^) = -Im^(g^/i^^ 

TT 



(An) 





Cln 


s 






A+ — 




+ x 


} 


I TT 











(2.11) 



with s = q^. It is straightforward to see that $(s, /i^) satisfies a homogeneous renormal- 
ization group equation (RGE), 



fi $(s, /i^) + 7 $ + $ 7^ 
djji 



0. 



(2.12) 



This implies that 

$(.) ^ c^{^?)^s,^^'')c\^^') (2.13) 

is a renormalization group invariant quantity, with C(yU^) being the Wilson-coefficient 
function of the 4-quark operators, which satisfies the RGE 

d 



V djji J 



(2.14) 



At the next-to-leading order, the coefficient function for AS" = 1 operators can be found 
in refs. The scale- and scheme-independence of $ should be clear, because this 

function is just proportional to the physical spectral function ;^Im\E'^'^^^(s). The scheme 
independence of Ti.^^^ is carried over to the 2-point function. 

We can easily sum up the next-to-leading logarithms in the 2-point function by setting 
/i^ = s, yielding 

o4 



e{s) 



{Any 



A + ^X 

TT 



and $(s) = C'{s)^{s)C*{s). (2.15) 



5 



Like the Wilson-coefficient function, also the spectral function at the next-to-leading 
order, in particular the matrices B and D, depend on the renormalization scheme. If the 
renormalization matrices in two schemes Za and Zf, are related by a finite shift. 



Za — Zh 



1 + — Ar 



(2.16) 



TT 

we find the following relation between X in the two schemes, 

Xb = Xa + ArA + AAr^ . (2.17) 

Together with the scheme dependence of the Wilson-coefficient functions (eq. (3.6) of 
ref. m), 



1 '^^ A T 

1 Ar 

TT 



a(/i'), (2.18) 



it is a trivial check that $(s) is indeed scheme-independent up to 0{a1). 

3 Current-current operators 

As an introductory example, we shall ffist calculate the 2-point function of the AS* = 1 
current-current operators, before embarking on the full set including penguins: 

Ql = iSaU/3)^^j^{ui3da)y^j^ , Q2 = isu)y_j^{ud)y_j^ , (3.1) 

where a, [3 denote colour indices (a, /5 = 1, . . . , A^) and the colour indices have been omitted 
for the colour singlet operator Q2. {V — A) refers to 7^(1 — 75). This basis closes under 
renormalization if penguin operators are neglected.! 

In the course of the calculation, it will become useful to also study 2-point functions 
of the Fierz-transformed operators 

Ql = (srf)v-A (^w)v-A ' = (^aC?/3)v-A (^/3Wa)v_A ' (3-2) 



and mixtures of the bases ( p.l| ) and ( p.2| ). Since these mixtures do not close under renor- 
malization, we will have to include evanescent operators. This will be discussed in detail 
below. 

The calculation of the 2-point function requires the evaluation of the leading order 
3-loop diagrams of fig. 1 and the next-to- leading 4-loop diagrams of fig. 2. The results 
of this evaluation are summarized in tables |] and H, and will be discussed in great detail 



^For the HV scheme, 7^ has to be taken in 4 dimensions. 
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in the following. A straightforward inspection reveals that the topology 2e contains traces 
with an odd number of 75's. Actually, there are two diagrams of type 2e: one with the 
fermion lines in the upper and lower loop circulating in opposite directions, denoted by 
2e, and one with the same direction, denoted by 2e'. These cannot directly be calculated 
in renormalization schemes with a naively anticommuting 75. For this reason, and also to 
make direct contact with the NLO calculation of the Wilson-coefficient function 
we shall perform the calculation with a non- anticommuting 75, originally due to 't Hooft 
and Veltman [^^|2^,0, and in addition we present a way to nevertheless obtain results in 
the NDR scheme. 



Table 1: Results for the lowest-order diagrams of fig. 1. 



Diagram 


la 


lb 


A 


4 

45 


4 

45 


j^NDR 


653 
675 


593 

675 


qHV 


1637 
1575 


1637 
1575 



Table 2: Results for the 0{as) diagrams of fig. 2 (Feynman gauge). 



Diagram 


2a 


2b 


2c 


2d 


2e 


2e' 


2f 


2g 


2g' 


2g" 


C 


1 

45 


1 

45 


2 

45 


2 

45 


8 

45 


2 

45 


8 

45 


4 

135 


4 

135 


4 

135 


jjNDR 


61 
180 


19 
60 


67 
90 


7 

10 


26 
9 


67 
90 


122 
45 


31 
81 


167 
405 


179 
405 


jjHV 


1349 
3780 


1349 
3780 


1643 
1890 


1643 
1890 


2866 
945 


1643 
1890 


2866 
945 


263 
567 


263 
567 


263 
567 



3.1 Current-current operators in the HV scheme 

Since the calculation is more transparent in the HV scheme, let us begin with this case. We 
shall denote with a matrix element of the general matrix eq. ( |2.2| ), corresponding 
to the 2-point function of the operators Qi and Qj. Then the three entries for the 2x2 
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matrix for Qi and Q2 (recall that \1/ is symmetric) are given by 



= iV2^ia + iV2c44^2a + 2^2e'] , (3.3) 

^12 = iV^la + ArC/[4^2a + 2^2c + 2^2e + 2^2e'] , (3.4) 
^22 = N^^,a + N''Cf[A^2a + 2^2c]+NCf<i/2g, (3.5) 

where Cf = [N'^ — 1)/2N. The contributions to the 2-point function from a given diagram, 
^diag, can be obtained by inserting into eq. ( |2.2| ) the relevant entries of tables [I| and 0. 
Including the corresponding colour factors and multiplicities which can be read off from 
eqs. (pD-d^D, we obtain the 2 x 2 matrices A, B^^ , C, and D^^: 

D«v |. (3.7) 

315 ^ \^ -1321 42iV ' ^ ' 

Inserting these results into eq. (|2.9|) , we arrive at 

X""- ^^NcJ''" -'^']. (3.8) 

225 ^ \^ -121 15iV J ^ ' 

Although at this stage a statement about the size of the radiative corrections is scheme- 
dependent, let us nevertheless perform this exercise. Taking as{s)/TT ^ 0.1, from eq. (|2.15D 




we find a moderate 20% correction in the diagonal, but the off-diagonal terms are almost a 
factor of 2 compared with the leading term. This already gives an indication of huge radia- 
tive corrections in the final result. However, note that these contributions are subleading 
in an expansion in 1/A^. 

Several technical remarks on the calculation so far are in order: 

i) Up to now, we only dealt with the current-current operators Qi and Q2- In this case, 
the penguin type contribution of diagram 2g in eq. (p.5|) has been omitted for consistency. 
It will be taken into account in the full result including penguin operators. 

ii) Naively, the HV scheme breaks some Ward-identities, e.g., the weak current is not 
conserved. We can enforce conservation of the weak current by performing a finite renor- 
malization which results in a shift for D^^ . This shift is given by —2CfA, and has been 
incorporated! into eq. ( |3.7| ). 



•^See also the discussion in refs. 
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iii) In the course of the calculation of the anomalous dimension matrix for 4-quark 
operators tensor structures with six 7-matrices appear which have to be projected 

onto the physical subset of operators. One example for the projections used in refs. 0,^] is 

7m7.7a(1-75)®7^7V(1-75) 4 (4 - 5) 7^(1 - 75) ® 7^(1 - 75) • (3.9) 

Generally speaking, at 0{6) this projection is arbitrary and a specification of the projection 
has to be added to the definition of the renormalization scheme. The projections used in 
refs. 1^, H have been chosen such that Fierz-relations in the current-current sector are 
preserved. This were not the case for an arbitrary projection. As a check, we have also 
performed the calculation with an arbitrary projection, and have verified that scheme- 
invariant quantities are indeed independent of this choice, as they should. 

Now, the 2-point functions have to be calculated in accord with the calculation of the 
anomalous dimensions. This means, we first have to calculate the radiative correction to 
either of the operators, then perform the projection onto the physical basis, and finally 
insert the resulting expression into the 2-point function. The only place where this treat- 
ment gives a different result, compared to a naive evaluation of the 2-point function (given 
the above choice of projection), is in diagrams 2e and 2f in the HV scheme. In the NDR 
scheme the naive calculation immediately yields the correct result for all diagrams except 
for the problems with 75 in 2e and 2e'. 

From tables |1| and it can be seen immediately that the result in the HV scheme 
respects Fierz-symmetry. Namely, the entries for the Fierz-conjugated diagrams (la, lb), 
(2a, 2b), (2c, 2d, 2e'), (2e, 2f), and (2g, 2g', 2g") are equal. In the case of the NDR scheme, 
we have the relations 

'^ZU,2f = a+eX''l2c,2e, (3.10) 

= vl/^^^^, (3.11) 
^^/^ = {l + e)^g?^ = {l + sY^^j:!^. (3.12) 

3.2 Current-current operators in the NDR scheme 

As was already mentioned above, diagrams 2e and 2e' contain traces with an odd number 
of 75's, and thus a direct evaluation in the NDR scheme is not possible. We can, however, 
use a "trick" in order to circumvent this problem. The trouble stems from the fact that 
Qi is in the colour non-singlet form. For 2-point functions of only colour singlet operators 
the problematic diagrams 2e and 2e' do not arise. Therefore, a solution lies in choosing 
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the basis (Qi, Q2), enabling one to calculate all diagrams without 75-problems |]26[. The 
price to pay is the fact that this basis is no longer closed under renormalization, and we 
have to explicitly add two evanescent operators Ei = Qi — Qi and E2 = Q2 — Q2- 
Now, the original and the new basis read 



Q 



Qi, Q2, El, E2) 



and 



Q = ( Q11 Q2, El, E2) . 



(3.13) 



The 1-loop renormalization of these two bases can be obtained by calculating matrix 
elements of Q and Q between free quark states. A straightforward computation leads to 



< > = MQ 



M 



Tree 



with 



/ -l/N 



1 + 



3 



V 



1 







1 

-l/N 







-1 

-1/A^ 
1 



1 \ 


1 

-l/N ) 



(3.14) 



The corresponding matrix M for the Q-basis is the same except for zeros in the entries 
(1,4) and (2,3). Thus, the finite shift Ar of eq. ( p.l7| ), mediating between the Q-basis and 
the Q-basis is given by 

/ 



Ar 










1 






-1 \ 










(3.15) 



For the 1-loop anomalous dimension matrices we find 
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(1) 



2Z\ 



(1) 



I -1/^ 
1 



V 



1 

-1/iV 







-1 

-l/N 
I 



I \ 


1 

-l/N ) 



(3.16) 



and 7*-^^ is the same with the entries (1,4) and (2,3) again being zero. 

We are now in a position to calculate all matrices in the NDR scheme. A and C are 
not scheme-dependent and therefore given by eqs. ( ^.61 ) and (|3.7D, with all other entries 
being zero. The matrices B^^^, B^^^, and D^^^, can be obtained from tables |^ and 0. 
The result is 



B 



NDR 



653 

"eo" 



A. 



D 



NDR 



45 ^ 



QN -175 
-175 6A^ 



(3.17) 
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B 



NDR 



675 



N 



f 653Ar 


593 


60 AT 


60 




593 


653iV 


-60 


-60iV 






-60 


120iV 


120 




V 60 


-60A^ 


120 


120A^ 


/ 



(3.18) 



where for B^^^ and D^^^ we only need the projection on the physical subspace. Next, 
we calculate the combination P = {•y^^^'B + By^^^'^)/2 which appears in eq. (|2T 



jNDR 



653 
225 



NCf 



1 

1 



yNDR 



593 
225 



NCf 



1 

1 



(3.19) 



In addition, the contribution from the shift Ar in eq. ( p.l7| ), ArA + AAr , vanishes. 
Combining everything, we obtain for D'^^^, 



D 



NDR 



D 



NDR 



P 



NDR 



P 



NDR 



— NCf 
45 ^ 



QN -187 
-187 6iV 



(3.20) 



Using this result together with eqs. ( p^ ) and 
D2gP^ of tab. H, which were not calculable directly. Finally, we have 



we can deduce the entries I^^'^^ and 



X 



NDR 



75 



NCf 



5N -47 
-47 5N 



(3.21) 



As a test for this result, we can use the relation ( p.l7| ) between the HV and the NDR 
scheme. The corresponding matrix Ar for this case can be obtained from eq. (3.9) of 
ref. § and eq. (4.25) of g: 



Ar 



A/N 1 
1 -1/A^ 



(3.22) 



It turns out that the relation ( |2.17| ) is indeed satisfied, providing a strong check of our 
result. 



3.3 The diagonal basis 



Further insight in our results for the current-current operators can be gained by trans- 
forming to the diagonal basis Q± = {Q2 ± (5i)/2 |1|,|31- Following ref. J^, we define the 
scheme-invariant operators 



Q± = 



l + ^B 



± 



(3.23) 
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where the scheme-dependent coefficients B± are given by 



B 



HV 



7 



1 



B 



NDR 
± 



The 2-point functions in that basis are exphcitly scheme-independent and read 



2B. 



TC 



^11 ±^12 



(3.24) 



(3.25) 



if we take advantage of the result \l/2c = ^2e'- 

Using this expression together with eq. ( |2.11| ), the spectral functions $±± 
turn out to be 

„4 



-Imvlf 



$±±(s,/i^) = e{s) 



(4vr) 









s 




Hi 






I 71 









3 101 43 1 
^4 ^^^lOiV 



(3.26) 



with A± = 2N{N ± l)/45. The coefficient of the logarithm is, of course, just equal to the 
leading-order anomalous dimensions of Q±, 'y±:\ 

The corresponding coefficient functions C±(/i^, M^) have been calculated in refs. |1|,0. 
We find it convenient to split up the Wilson-coefficients in two factors, one solely depending 
on and the other on M^: C±{^'^, ^w) ~ C'±(/i^) C±{M'^). The two factors are given 
by 



Air 



(3.27) 



The NLO correction R± can be found in ref. 0. 

Using this result, we are in a position to form the scale- independent spectral functions 
of eq. (gJEi), <l±±(s) = Cl{M^)Cl{s)^i±{sy. 



a,(s)2^±'/'3iC|(M, 



(47r)f 



w) 



71 



Setting / = 3, we find for the NLO contributions X±: 



X4 



121 

540 
121 



30917 
16200 



,,4 22997 . 



11173 
5400 

8143 



781 
6OO' 

641 



5 

12 



540 16200 
For = 3 the two spectral functions simplify to 



-A^^ + 

5400 600 12 



2 ^ 
W) 



3649 aAs) 



2 -1 

W) 



1620 77 



9139 as{s) 



(3.28) 

(3.29) 
(3.30) 

(3.31) 
(3.32) 
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Let us comment briefly on the implications of our results. Again taking a;s(s)/7r ^ 0.1, 

at the NLO we find a moderate suppression of by roughly 20%, whereas $ acquires 

a huge enhancement on the order of 100%, including the coefficient functions at Mw, 
C±(M^), which only have a minor effect. Because $++ solely receives contributions from 

AI = 3/2, and $ is a mixture of both A/ = 1/2 and A/ = 3/2, this pattern of the 

radiative corrections entails a strong enhancement of the AI = 1/2 amplitude. Hence, we 
are provided with a very promising picture for the emergence of the A J = 1/2-rule. 

Analyzing our result from the point of view of the large- expansion [l^,|T^, we see 

that at leading-order, the corrections to $++ and $ are equal [^4|,^, meaning that 

the AI = 1/2-rule is missed completely. This situation is partly remedied if the large 
next-to-leading corrections in 1/A^ are taken into account. 



4 Full result including penguins 

To obtain the complete result, we have to add to our basis the following four penguin 
operators 



^"/SJv-aXI 
Q 



(Sad, 



aJy-A ' 



(4.1) 



which arise from the current-current operators in the process of renormalization. The 
corresponding Fierz-transformed operators, again being needed for the calculation in the 
NDR scheme can be found in ref. p|. 

The 2-point functions for {V — A) ® {V — A) operators including Qs and Q4 are given 

by 



^^13 

^^23 
^33 
^44 



4^26 + 2^ 



2d 



N'^ib + NCf 4^26 + 4*2d + 2^ 



'2f 



^U + 2NCf^2g, ^24 = ^13 + /A^C/^2c;' , 

/^ll + 2^23, ^34 = /^12 + 2^'24, 

f^n + 2^14 + fNCf^2g" . 



(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 



In the expression for \E'33, we have used the relation \E'2c = ^2e'- 

The 2-point functions with insertions of the (V — A) (V + A) operators Q5 and Qq 
can be calculated analogously. For simplicity, we don't give their contributions in detail. 
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but just state the final results. In the HV scheme, we find: 



A 



B 



HV 



N 







( N 


1 


1 


n n \ 








1 


N 1 


A^ 







A 

4 






1 /A^ + 2 2A^ + / 


n n 
U U 




45 




1 


N 2N + f fN + 2 




















fN f 








I 








f fN J 




1637Af 


1637 


1637Ar 


1637 


-2N 


-2 


1575 


1575 


1575 


1575 


315 


315 


1637 


1637Ar 


1637 


1637Af 


-2 


-2N 


1575 


1575 


1575 


1575 


315 


315 


1637Af 


1637 


1637{/Af+2) 


1637{2Af+/) 


-2{/Af+2) 


-2{2Af+/) 


1575 


1575 


1575 


1575 


315 


315 


1637 


1637Ar 


1637(2Af+/) 


1637(/Af+2) 


-2{2N+f) 


-2(/Af+2) 


1575 


1575 


1575 


1575 


315 


315 


-2N 


-2 




-2(/7V+2) 


-2{2N+f) 


1637/Af 4 


1637/ 4Af 


315 


315 


315 


315 


1575 315 


1575 315 


-2 


-2N 


-2{2Af+/) 


-2{fN+2) 


1637/ 4:N 


1637/Af 4 


315 


315 


315 


315 


1575 315 


1575 315 



(4.7) 



(4.8) 

The matrix C satisfies the relation (p.lO|) with 7^^^ given in ref. p|, and we don't list 
it explicitly. The matrix D^^ has been relegated to the appendix. Inserting these into 
, we obtain for X^^: 



eq. 



/ 27V 
15 



X 



HV 



NCf 



-242 



225 



-242 
225 

2N 
15 

-242 
225 





2N 
15 



242 
2025 



2N 
15 

-2662 
2025 



-242 
225 

2N _ 242/ 
15 2025 



-2662 



2025 



2N 
15 



242/ 
2025 



2fN_ 

15 

47V _ 
15 



_ 5324 
2025 

2662/ 
2025 



47V 
15 



2fN 
15 



V 





-242/ 
2025 





-484/ 
2025 



_ 2662/ 
2025 

242/^ _ 
2025 



-242/2 



484 
225 








2/7V 
15 

38/ 
25 





-242/ 
2025 

-484/ 
2025 

-242/2 
2025 

38/ 
25 



322/7V 
225 



242/2 
2025 



2025 25 225 2025 / 

(4.9) 

The treatment to work around the 75 problem in the NDR scheme for the full basis 
parallels the method used in the current-current case. We can choose the basis used in 
ref. ||2^, (Qi, Q2, Q3, Qa, Q5, Qe), which does not contain colour non-singlet operators, 
thus not posing 75 problems for the computation of the 2-point function. However, because 
this basis does not close under renormalization, this time, for each of the six operators we 
have to add to the basis an evanescent operator, implying that at intermediate steps of the 
calculation we have to work with 12 x 12 matrices. We skip the unilluminating details of 
this computation and just present our results. 
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As already stated, the matrices A and C do not depend on the scheme, and agree with 
the HV case. For B^^^ we find 



B 



NDR 



N 



V 



653Af 


653 


593Ar 


593 


675 


675 


675 


675 


653 


653Af 


593 


593Af 


675 


675 


675 


675 


593Ar 


593 


653/Af+1186 


1186Af+653/ 


675 


675 


675 


675 


593 


593Af 


1186Af+653/ 


653/Af+1186 


675 


675 


675 


675 



















































653/Af 


653/ 


675 


675 


653/ 


653/Af 


675 


675 



(4.10) 



and D^^^ can again be found in the appendix. From these results we deduce for X^^^: 



X 



NDR 



NCf 



/ 


2N 


15 




-94 




75 




2N 




15 




-94 




75 







\ 






-94 
75 



2N 
15 



182 
2025 



2JV 
15 

-2902 
2025 



-2902 



2025 



2N 
15 



212/ 
2025 



2fN_ 
15 

4N _ 

15 



_ 5804 
2025 

2962/ 
2025 



2N 
15 

4N 
15 



2fN_ 

15 





-212/ 
2025 





-424/ 
2025 



-94 
75 

_ 212/ 
2025 

_ 2962/ 
2025 

242/^ 
i) 

-242/2 



188 
75 



2025 









2fN_ 

15 

74/ 
75 





-212/ 
2025 

-424/ 
2025 

-242/2 
2025 

74/ 
75 



94/Ar 
75 



242/2 



2025 



(4.11) 

The expressions for X^^ and X^^^ given above again do satisfy the relation ( p.l7|) 
with Ar taken from ref. IHl a 



Ar 



eq. 


(4.25) of i 
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6N 
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6N 
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6 








2 
3N 


2 

3 


1 

3N 


1 
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N 




















3 
N 


3 














2 


N - — 



\ 



(4.12) 



V 

This test of our results provides us with great confidence as to their correctness. 
Let us make a few observations on the results thus obtained: 

• Comparing the NLO matrices X^^^ and X^^ to the LO matrix A, we find huge 
corrections on the order of 200% for the entries (1,2), (1,4), (2,3), (5,6), and (6,6). 
The difference between the NDR and HV scheme is small with respect to the large 
absolute value of the corrections. All other entries have moderate corrections ^ 50%. 
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If the number of flavours / = 3, in the HV scheme we have the operator relation 
Qi = Q^i + Q2 — Qi- This leads to the following relation for the two point functions: 
'^^^ = "iff^j^ + '^2i^ — "^li^ . The relation is satisfied by our result, providing an 
additional test. In the NDR scheme, the relation is broken by 0{af) corrections.@ 

Due to the factorization property of the Qq operator in the large- limit, in this limit 
\E'66 can be related to a convolution of two 2-point functions for scalar currents p6 |. 
This relation holds for our result but we shall return to this point in sect. 7. 



5 Numerical Results 

In this section, we shall provide the reader with a brief discussion of the numerical implica- 
tions of our results, and postpone a thorough phenomenological analysis to a forthcoming 
publication. 

Following the notation of refs. the Wilson-coefficient functions for AS* = 1 weak 

processes can be decomposed as C{s) = z{s) + Ty{s), where r = — (Vj^V^^) / (y^^V*^). 
The coefficient function z{s) governs the real part of the effective Hamiltonian, and y{s), 
parametrizes the imaginary part and governs e.g. the measure for direct CP-violation in 
the i^'-system, e'/e. We thus have two different quantities with the help of which we can 
form the scale- and scheme-invariant combination of eq. ( p.l5|) . Let us denote these 
two functions by: 

= z''{s)^s)z{s); %{s) = y^{s)^s)y{s). (5.1) 
For the numerical analysis, we shall consider the range Q = y/s = 1—3 GeV, appearing 



as a natural scale for a QCD sum rule analysis of the i^-system |25]. In this range, the 
coefficient functions have the following structure: 

• Above the charm threshold rric, z^ — zq vanish, and $2 is only given as a product of 
{zi, Z2) and the current- current part of Below rric, penguins are generated from 
the operator mixing, but the coefficient z^ — Zq still remain small above 1 GeV, such 
that in the whole range is dominated by current-current operators. 

• In the case of y{s), only y?., — y% are non- vanishing in the whole range considered. The 
coefficient for the Qq operator y^ dominates, but the other penguin operators also 
give noticeable contributions. 



*See also sect. 4.5 of rcf. 
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Since in this work we are mainly interested in the size of the radiative corrections to 
the effective Hamiltonian, we write as 

^^,yis) = (5.2) 

where the superscripts (0) and (1) refer to the leading as well as next-to-leading order 
respectively. In fig. 3, we plot the ratios and $W for Aj^^ = 200, 300, and 

400 MeV. For simplicity, in fig. 3, we have not included the quark threshold at m^, but we 
work in a theory with / = 3 up to 3 GeV. We have checked that these threshold effects 
only cause a small change in the NLO correction. Of course, as expected, the values for 
^z,y in the HV and NDR scheme exactly agree. 

From fig. 3, we can see that in the region Q = 1 — 3 GeV, and for a central value 
^MS — 300 MeV, the radiative QCD correction to $2 ranges approximately between 40% 
and 120%, whereas in the case of $y we find a correction on the order of 100%-240%. 
Because the 2-point function is constructed as the square of the effective Hamiltonian, the 
actual corrections to Ti.^^^ are only about half the corrections to the 2-point function. 
Therefore, the perturbative QCD correction to the real part of the effective Hamiltonian 
turns out to be 20%-60%, and for the imaginary part 50%-120%. 

In sect. 3.3, we have demonstrated that the large corrections correspond to the 
A/ = 1/2 part of the effective weak Hamiltonian!. The corrections to the A/ = 3/2 
part are identical to the ones in the AS* = 2 correlator (both operators are in the same 
representation of the chiral group), which, as shown in sect. 3.3 and the next section, are 
quite moderate and negative. This implies that for the A/ = 1/2-rule in K it tc decays, 
we receive an additional large and positive contribution, bringing theoretical calculations 
closer to the experimental value. 

The calculation of the imaginary part of Ti-es, no longer retains perturbative character, 
because of the large corrections. Nevertheless, this does not completely spoil existing 
determinations of weak matrix elements in the framework of the 1/A^ expansion or chiral 
perturbation theory, for there, to a given order in 1/A^ or the chiral expansion, the largest 
corrections in ag are completely summed to all orders. 

^This has also been extensively discussed in ref. |2q]. 
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6 The A5 = 2 operator 



For the case of AS = 2 transitions, things are somewhat simpler because there is only one 
operator. We take this operator to be 

QaS=2 = ^ [ (srf)v_A (sC?)v-A + iSadp)y_^ (%rfa)v_A ] • (^-l) 

This definition might seem unfamiliar, but we shall discuss in the following, why it is 
convenient. First, let us note that it renormalizes into itself, even in a general dimension!, 
and it is obviously Fierz-symmetric. Apart from the quark content it has the same structure 
as Q+. 

Collecting the contributing diagrams, and making use of the relation the 2-point 

function of Qa5=2 is given by 

^AS=2 = ^11 + ^12 + ^13 + ^14- (6.2) 

From this expression, we obtain the following values for A, B, and X: 

Aas.^ = ^N{N+i), (6,3) 

BiL, = - 1^ iv(iV + l) ; B^ia = - ^ iv(iV + l) , (6.4) 

As was already remarked at the end of the last section, up to a multiplicity factor 4, 
these quantities agree with the corresponding expressions for \E'++ of sect. 3.3, if we would 
refrain from performing the rotation of eq. ( ^.23| ) to a scheme-invariant basis. Cas=2 again 
respects eq. ( p.lO|) , and the relation between schemes ( p. 17] ) is also satisfied with 



ArAS=2 = o ( 1 - ^ ) , (6.6) 



being easily obtained from ref. f^. 

In the HV scheme, we could have worked with the operator 

Oas=2 = (sc?)v_A (sc?)v-A ' (6-7) 

in 4 dimensions being equivalent to Qas=2, since in HV Fierz-symmetry is respected for 
current-current operators, and Oa5=2 renormalizes into itself. This is not true for the NDR 



^ Apart from evanescent terms which are taken care of by the projection discussed in sect. 3.1 



18 



scheme and we would have to go through a similar procedure as described in sect. 3.2. 
Namely, augmenting the basis by an evanescent operator Eas=2 = Oas=2 — 0^3=2 and 
including this operator for renormalization, because it induces additional contributions to 
the physical subspace. We have checked that this treatment leads to the same value for 
-^As=2- Let us point out that this only concerns the quantities B^g3^ and D^g32. The 
anomalous dimensions of Qas=2 and Oas=2 in the NDR scheme agree even at NLO. 

In order to be able to form the scheme- independent combination $a5=2(s), we need 
the Wilson-coefficient function for AS" = 2 processes at NLO. It can be obtained from 
refs. [^1^ for internal top and charm quark exchange in the box-diagram. The mixed 
charm-top contribution is not yet available at NLO. However, here we shall not pursue 
this any further, but use the strategy of ref. of defining a scale- and scheme-invariant 
operator for AS* = 2. This operator is given by 



Q 



AS=2 



47r 



Q 



AS=2 



where 7^5=2 — 7+"* is the LO anomalous dimension of Qas=2- The finite NLO correction 
Z depends on the scheme, and can be found in [^. The matrix element of this operator 
is directly parametrized in terms of the scheme-invariant S-parameter for — K'^- 
mixing. 

Calculating the spectral function $a5=2('5) for Q/\s=2, we obtain 

„4 



'^AS=2is) 



(An) 



^AS=2 H ^AS=2 

IT 



(6.9) 



with Xa5=2 = 4X_|_, and X+ being given in eq. (|3.29|) . This function is explicitly scheme- 
invariant, and for = 3 takes the form 



AS=2[S) 



9{s] 



32 



15 (47r) 



a.As 



-4/9 



3649 aJs) 



1620 TT 



(6.10) 



Because both, and Qas=2, have the same chiral representation, as expected, apart from 
a global factor, their spectral functions agree. We observe that the NLO QCD-correction 
is negative and on the order of 20%, for as{s)/Tr ^ 0.1. 



7 Discussion 



Our work improves and completes the 2-point function evaluation of ref. [^] with two ma- 
jor additions: the recently calculated NLO corrections to the Wilson-coefficient functions 
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have been taken into account and, moreover, we have incorporated the missing contribu- 
tions from evanescent operators. The final results are then renormalization scheme- and 
scale- independent at the NLO, and, therefore, constitute the first complete calculation of 
weak non-leptonic observables at the NLO, without any hadronic ambiguity. 



It is worthwhile to make a comparison with the results of ref. p6|. In this work, the 
NDR scheme was used and the calculation of the 2-point function for the AS = 1 case 
was performed in the basis 

Q = Q2, Q3, Q4, Q5, 4), (7.1) 

not being plagued by problems with 75. The topologies present in that calculation were 
la, lb, 2a, 2b, 2c, 2d, 2f, and 2g. We have reproduced the results for all these diagrams 
and we fully agree with ref. [^]. However, as already remarked in sects. 3 and 4, in the 
NDR scheme the basis (17. 11) does not close under renormalization, and we had to add 



contributions from evanescent operators which were not included in ||2^. In the notation 
of sect. 3, these evanescent contributions shift the matrices B^^^ and D^^^ to B^^^ and 
jjNDR^ and therefore the final correction X^^^ gets changed. In addition, the relation 



Q4 = Qs + Q2 — Qi for / = 3 was used in ref. [g6| to eliminate Q4. As it stands this 
relation is valid on the operator level. 

Removing the entries for from our matrices B^^^, D^^^, and X^^^, and setting 
/ = 3, one can easily see the differences with the results of ref. |2^. The evanescent 
contributions have changed the entries (1,2), (1,3), and (6,6) of B^^^, and (1,2) and (6,6) 
in D^^^. The differences in B^^^ propagate via the matrix products of eq. ( p.9|) into 
most entries of X^^^: all entries except for (5,5) and the trivial zeros in (1,5), (2,5), and 
(3,5) are different. The final numerical differences are however not big, since the most 
sizeable corrections were already included in the original calculation of ref. |2^ . 

In the large- limit the operator Qg factorizes in the product of two current operators. 
Therefore, the 2-point function ^fge can be calculated as a convolution of two current- 
correlators.0 This was used in ||26| as a check for ^fge at the leading order in 1/A^. In 



fact, the large- limit result of ref. was already a full NLO calculation, since the 
corresponding anomalous dimension was already known at the NLO (it is related to 
the quark-mass anomalous dimension in the large- A^ limit) and it was correctly taken into 
account. Although our results for B^^^^ and D^^^^ show a discrepancy to the result 



of ||26| even at the leading order in 1/A^, the combination Xqq only deviates from the 



result of |2^ by sub leading terms in the 1/A^-expansion, hence fulfilling the test in both 



^For details see ref. 
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cases. The differences at intermediate steps of the calculation {Bqq^^ and D^/^^) stem 
from the fact that a different form of the operator {Qq or Qq) is being used, but the final 
physical result is of course identical. 

For the AS* = 2 operator the evanescent contributions result in changes in all next-to- 
leading quantities -Ba^S^ ^^3=2^ ^"^^ ^^.3=2- The coefficient of the NLO correction to 
^/^s=2, eq. (|6.10|) , was found to be —1217/810 in ref. [Q, compared to —3649/1620 in our 
case. This correction was used in ref. I^S] for a sum rule determination of B^- The effect 
of our new result would be to slightly further reduce the value of Bk obtained in ref. p5| . 



Qualitatively, the conclusions of ref. ^6] remain unchanged. At the NLO, the AJ = 1/2 
piece of the AS* = 1 effective Hamiltonian gets a huge positive correction, while the gluonic 
effects in the A/ = 3/2 (and AS* = 2) operator are moderate and negative. Together with 
the previously known enhancement of the Wilson-coefficient this provides a very 

suggestive explanation of the observed enhancement of A/ = 1/2 transitions in K decays. 

As explicitly shown in fig. 3, the NLO gluonic contributions are even more important 
in the CP-violating piece of the weak AS* = 1 Hamiltonian; the reason being that this 
part is dominated by the penguin operator Qe, which gets the largest correction. As 
shown in ref. [^, this enhancement is further reinforcedi at 0{a1)^ indicating a blow-up 
of the perturbative series in this case. Fortunately, this non-perturbative character does 
not completely prevent the feasibility of a reliable determination of CP-violating effects, 
since these leading contributions in 1/A^ can be resummed to all orders. 

In a series of articles a more phenomenological description of the A/ = 

1/2-rule was advocated. The key idea is to rewrite the 4-quark operators as a product of 
diquark-anti-diquark operators by means of Fierz-transformations, and treating the diquark 
as an effective particle, similar to the constituent quarks. The important observation then 
lies in the dominance of pseudoscalar diquark matrix elements for low momentum transfer 
over axialvector meson matrix elements which are proportional to the momentum. All 
low energy decays in which diquarks can participate show the enhancement of A/ = 1/2 
amplitudes and a surprisingly good description of those processes was obtained. 

An inspection of our results at the diagrammatic level reveals the following pattern for 
the origin of large corrections: all 2-point functions (except those which vanish at lowest 
order) receive contributions from the self-energy diagrams 2a or 2b, as well as from the 
quark-antiquark vertex corrections 2c, 2d, or 2e'. These contributions cancel to a fair 



^Thanks to the factorization property of the penguin operator in the large- iV hmit, the 0{a1) correction 
is also known in this limit. 
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amount: 



C2-2a+2c — C2-2a+2e' — C2-2b+2d — , (7-2) 

j-^NDR _ p,NDR _ p,NDR _ ^ ^7 o\ 

-^2-2a+2c ~ -'-^2-2a+2e' ~ ^2-2b+2d ~ ' \'-'^) 

r,HV _ T^HV _ p,HV _ ]_ 

-^2-2a+2c ~ -^^2-20+26' ~ ^2-2b+2d ~ 45 ' \' -^J 

The quark-quark and antiquark-antiquark correlation diagrams 2e or 2f already by them- 
selves are the biggest terms, and due to the partial cancellation of self-energy and current- 
vertex diagrams, we receive large corrections wherever quark-quark correlations can con- 
tribute. Note that these diagrams are subleading in The penguin diagrams 2g, 2g', 
and 2g" generally have small impact on the 2-point functions. 

This structure of the radiative corrections to 2-point functions of AS = 1 and AS* = 2 
operators allows for a deeper understanding, why the description of non-leptonic weak 
decays in terms of diquarks was so successful as far as the A/ = 1/2-rule is concerned. 
In this framework, by working with effective diquarks, the quark-quark correlations were 
phenomenologically summed up to all orders in the strong coupling. Since these are the 
dominant corrections to the 2-point functions, summing them up provides us with a very 
good physical picture of the underlying QCD dynamics. As such, though, the statement in 
question, like the diquark-current itself, i is gauge-dependent. In fact, the gauge-invariant 
combinations of diagrams are (2 ■ 2a + 2c) , (2 ■ 2a + 2e') , (2e + 2e') , as well as their cor- 
responding Fierz-conjugates. However, now the gauge- independent combination involving 
the quark-quark correlations dominates the other terms even more drastically (by one order 
of magnitude): 

2 

^^26+26' = C2d+2f = ' {^■'^) 



T^NDR _ j^NDR _ p\ 

^2e+2e' " , ^2d+2/ " ' ^'■'^> 



r)Hv _ p)Hv _ 1363 

-^2e+2e' ~ ^2d+2f ~ ggg • K' ■' J 

A full QCD calculation has been possible because of the inclusive character of the 
defined 2-point functions. Although only qualitative conclusions can be directly extracted 
from these results, they are certainly important since they rigorously point to the QCD 
origin of the infamous A/ = 1/2-rule, and, moreover, provide valuable information on the 



'See the related discussion in ref. [R5| 
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relative importance of the different operators, which can be very helpful to attempt more 
pragmatic calculations. Obviously, a direct application of our results would be the use of 
dispersion relations to extract "more exclusive" information from the 2-point functions, 
following the methods developed in refs. |20,|21],E2|,p3lEl,l2^,|26|. We plan to investigate 



this and other possible phenomenological applications in the future. 
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Figure 1: Leading order diagrams for the 2-point function. 
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